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(57) Abstract 

A control design as shown (figure 3C) having a set 
of linearized differential equations which are computed 
based on data from identification (340 of figure 3A). At 
step (351), a high order model is developed, at step (355), 
model reduction is performed by computing an approxi- 
mate model of the higher order model, of step (351), with 
a fewer number of states than the high order model which 
may provide an estimate of errors due to the approxima- 
tion. At step (356), estimator controller or controller syn- 
thesis may be performed including computing an estimator 
and/or controller that optimizes certain properties such as 
norm of the closed loop system. At step (357), the high 
order model is evaluated by computing multivariable fre- 
quency responses, step responses, responses to stochastic 
noise processes, eigenvalues, and/or measures of robust- 
ness for the closed loop system consisting of the high 
order model and the synthesized estimator controlier or 
controller. 
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ADAPTATION TO UNMEASURED VARIABLES 

REFERENCE TO RELATED APPLICATION 

This application claims the benefit of U.S. Provisional AppUcation No. 
60/115,165. 

FIELD OF THE INVENTION 

This invention relates to the field of control systems and, more specifically, to 
management of uncertainty in control systems that require adaptation. 

BACKGROUND 

The design of new systems is becoming increasingly complex because of 
reliance on high performance features requiring tight tolerances and more 
robustness to system variations, for example, in resonance and anti-resonance. 
Resonance is a phenomenon that occurs in different types of systems, including 
mechanical, electrical, and acoustic systems. Many industries, for example, the disk 
drive industry, rely on system component redeii^ in order to handle the 
resonances in their systems because current disk drive servo systems are unstable in 
the presence of resonances. 

In a disk drive, a servo system is typically used to position a head assembly 
relative to a disk in order to reproducibly read and write data onto a specific 
location of a disk. The industry trend is toward increasing data storage densities 
and decreasing data access times by packing data locations as close to one another 
as possible (i.e., increasing track density) and by faster positioning of head 
assemblies over the data locations (i.e., decreasing seek time). This, however, leads 
to problems with resonance that develops in servo systems having tight alignment 
tolerances. Mechanical resonances in servo systems are a result of complex 
vibration interactions in the disk drives. Whereas, electrical resonances and anti- 
resonances in channel transfer hmctions may be the result of transmission line ta^ 
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For example, in a magnetic disk drive, a head is flown over the disk surface 
on a thin air bearing. Mechanical vibrations in the drive due to air hirbulence may 
result in tracking (Le., head positioning) errors, thereby limiting achievable storage 
densities or resulting in increased access times to allow for resonance to dampen 
out of the system. 

Unmeasured variables may affect the operation of a drive and can cause 
variations in performance from unit-to-unit, as well as within a unit, during use of 
the drive. These variables may include, for examples: temperature, spingle bearing 
asperity, disk flutter induced non-repeatable run-out (NRRO), and variable 
repeatable run-out (RRO). 

One solution for reducing track errors is to redesign disk drive components 
that contribute to system resonance. Resonance may be engineered out by, for 
examples: designing smaller head assembUes, redesigning drive enclosures and 
mounts, and using disk materials that reduce disk flutter. Mechanically 
engineering out resonances, however, may be expensive and may delay getting a 
product to market. 

Another solution is to design a servo control system so that the bandwidth of 
the system is increased and the resonant frequency is close to or inside the 
bandwidth of the system. The system bandwidth can be generally described as the 
range of frequencies within which a system will respond satisfactorily with 
acceptable gain and phase response. The gain margin and the phai^e response are 
measures of the relative stability of a system. The phase response of a system 
changes rapidly near resonance. As such, increasing the system's bandwidth and 
moving the resonant frequency nearer to a system's cross-over frequency makes the 
servo system more sensitive to the gain and phase response of the resonance. 

One problem with these prior art servo systems is that as the systems 
bandwidth is increased, the performance of the system degrades due to unit-to-unii 
manufacturing variations as well as variations within a particular unit operated 
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under different conditions. The performance of a system may be defined, for 
example, in terms of root mean square (rms) in servo tracking errors as a response 
to given vibrations, measurement noise, and run-outs. As such, a resulting high 
bandwidth control system will have less robustness (i.e., less gain and phase 
margin). 

Another problem with prior solutions is that a large uncertainty results when 
attempting to design only a fixed disk drive confa-ol system using only one model 
that accounts for all system variations. A single conta-oUer that has to handle all the 
variations without adapting reqtiires high robusfaiess to large uncertainty in the 
control design models. Large model imcertainty leads to poorer achievable 
performance because the design has to achieve that performance over a much larger 
set of models- Therefore, a single contaroUer has to sacrifice performance to achieve 
robustness. 

As an example, it is difficult to predict structural frequencies of disk drive 
components using simulation models as unit-to-unit variations dominate variations 
in responses. In addition, the parameter estimation schemes used in prior art 
solutions have difficulty in accurately determining resonant modes in the presence 
of unmeasured stochastic and periodic disturbances with limited closed loop data. 
As such, the phase variations in the resulting transfer fimction of the drive control 
system are large, requiring high phase and gain margins nearing resonance. 
Furthermore, in-use variations such as those due to temperature and external 
vibrations create even more uncertainty. A larger uncertainty in the models used in 
conti-ol system design results in a lower performance in actual operation of the disk 
drive. 

As such, prior art solutions that do not adapt to the changes caused by 
unmeasured variables are limited by the accuracy of one control model of a system, 
and do not allow for increased contirol performance based on resource utilization^ 
design and manufacturing time. Altiaough, conventional parameter adaptive 
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schemes, for examples: Model Reference Adaptive Control (MRAC), Self-Tuning 
Regulator (SIR), and Pole-Placement, may be effective in low bandwidth, low 
systems, these methods have significant in-use robustness problems in high 
bandwidth systems (e.g., disk drive servo control systems). Adaptation to changes 
caused by unmeasured variables in resonant systems with such prior art methods is 
non-robust because prior art adaptive methods: 1) do not account for uncertainty in 
models; 2) attempt to handle adaptation to large in-use and unit-to-unit variations 
ignoring model uncertainty; and 3) do not attempt to match compute resource 
utilization with information that is available at design and manufacturing time. 
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SUMMARY OF THE INVENTION 

The present invention pertains to a method of generating a conta-ol system. 
The method may include generating a set of conta-ol models of the plant based on an 
operation of a prototype at a design time. The method may also include 
characterizing the operation of the plant at a manufarturing time and selecting a set 
of model contaroUers based on the characterizing of the operation of the plant at 
manufactiiring time. Each of tiie model conti-oUers corresponding to one of the 
control models. In one embodiment the method may also include selecting one of 
the model contiroUers at a run time based on an operating parameter. 

In one embodiment, the present invention includes a plant controller having 
model conti-oUers. Each of the model conta-oUers may be coupled to receive 
monitored signals and generate a control signal. The plant contiroUer may also 
include a data evaluator coupled to receive the monitored signals and generate an 
output. The data evaluator determines when sufficient information from the 
monitored signals is received. The data evaluator may be coupled to select the 
control signal of one of the model confa-ollers based on the output of the data 
evaluator. 

Additional featiires and advantages of the present invention will be apparent 
from the accompanying drawings and from the detailed description that follows. 
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BRIEF DESCRIFnON OF IHE DRAWINGS 

The present invention is iUustrated by way of example, and not by way of 
limitation, in the figures of the accompanying drawings and in which: 

Figure 1 illustrates an embodiment of a controller used with a disk drive 
system. 

Figure 2 illustrates an embodiment of the development of a control system. 
Figure 3A illustrates an embodiment a design cycle. 
Figure 3B illustrates an embodiment of a closed loop system on which data 
acquisition is performed. 

Figure 3C illustrates an embodiment of control design and control 
implementation. 

Figure 4 iUustrates an embodiment of a control model structure for system 
identification performed at design time. 

Figure 5 illustrates an embodiment of an imcertainty model. 

Figure 6 illustrates an embodiment of a multiple model controller. 

Figure 7 illustrates an embodiment of a hypothesis tester and test confidence 
indicator. 

Figure 8 illustrates an embodiment of a test confidence indicator. 
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DETAILED DESCRIFnON 

In the following description, numerous specific details are set forth such as 
examples of specific systems, components, parameters, etc. in order to provide a 
thorough understanding of the present invention. It will be obvious, however, to 
one skilled in the art that these specific details need not be employed to practice the 
present invention. In other instances, well-known materials or methods have not 
been described in detail in order to avoid unnecessarily obscuring .'he present 
invention. 

A method for adaptation to unmeasured variables in the design, 
manufachure, and contirol of a plant is described. A plant is a controlled element in 
a system that may be, for examples, a process, an article, and/or a machine. A 
conti-oller is employed to control particular values or conditions of the plant in 
response to factors that influence the behavior of the plant. In one embodiment, the 
contirol scheme includes: generating a group of control models of a plant at a 
design time based on characterizing and evaluating prototypes of tiie plant under 
various conditions; characterizing and evaluating tire-performance of the plant at a 
manufacturing time, and; selecting a sub-set of the group of control models based 
on the characterization of the plant at manufactixring time to control an individual 
plant during run time. In one embodiment, a metric of uncertainty is maintained 
that is used in selecting appropriate control models with specific robustaiess- 
performance tradeoffs. 

The contirol scheme described herein may be implemented with a wide 
variety of systems including acoustical, mechanical, and electrical systems. For 
example, the control scheme described herein may be used to provide control of 
resonance in a disk drive system as discussed below. It should be noted, however, 
that the control scheme is described in relation to a disk drive system only for 



-7- 



BNSDOCID: <WO_ .0041042A1 I > 



wo 00/41042 PCr/USOO/00387 



illustrative purposes and is not meant to be limited to control of resonance in disk 
drive systems. 

Figure 1 illustrates one embodiment of a controUer used with a disk drive 
system. The disk drive 100 includes a disk 110 having concentric tracks containing 
stored data, a head 120 for reading and writing the data on disk 110, an actuator 130 
to position the head over a desired track during read and write operations, and a 
servo controller 150 to control actuator 130. The movement of head .120 to a desired 
track is referred to as "seeking" and the maintenance of head 120 over the desired 
track is referred to as "settling." In one embodiment, achiator 130 is a voice coil 
motor (VCM). In other embodunents, actuator 130 may be other types of motors, 
for example, a stepper motor. 

A VCM contains a coU that is movable in a magnetic field of a magnetic 
stator. Applying current to the VCM causes the coil to move which, in him, causes 
head 120 to move in a radial direction across disk 110. The acceleration of the coil is 
proportional to the applied current. The track information is incorporated into disk 
110 and read by a servo control system using head 120. The servo controller 150 
uses the track information to maintain head 120 precisely aUgned over a track on 
disk 110 based on track setpoint signal 124. The track information is demodulated 
to generate a position error signal 121 indicating the position of head 120 away 
from the nearest track centerline. A position sensor signal (Ypos) is generated by 
adding the track position signal 123 to position error signal 121. The position 
sensor signal is used in the servo control system feedback loop to generate a control 
signal 155 to the VCM to move head 120 back to the track centerline. 

During operation of the VCM, resonant modes in the system cause variations 
in the coil current in the VCM, resulting in movement of head 120 away from the 
track centerline. Having resonant frequencies at or near the bandwidth of the 
system can lead to instabilities in the system. During seeking, resonances may not, 
have a significant affect on head 120 positioning as the primary contribution to the 
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total forces on head 120 comes from the current that is appUed to the actuator. 
However, during settling, resonance may become a primary contributor to head 
position error and must be quickly compensated for by the servo control system. 
The servo control system produces a current in the VCM in order to maintain head 
120 over the centerline of the track. In one embodiment, the servo controller 150 
includes a multiple model controller that enables the servo controller to adapt to a 
large range of resonance. Such a multiple model controller has the ability to handle 
a larger range of resonances than may be handled by prior art parameter adaptive 
controllers. 

Figure 2 illustrates an embodiment of the development of a control system. 
In one embodiment, the uncertainty in models characterizing a plant is successively 
reduced from design time 220 to manufacturing time 230, and then further reduced 
at run time 240, for the specific plant under specific operating conditions. At design 
time 240, a limited number of representative prototype units are available and are 
characterized at specific wide ranging operating conditions. In one embodiment, 
specific operating conditions may be, for example, specific values of temperature, 
humidity, age, and/or vibration spectrum, for a specific head, on a specific region 
of the platter. The available prototype drives are characterized at each of the 
specific conditions, using system identification techniques in a design cycle as 
discussed below in relation to Figure 3. Using uncertainty calculations for models 
and trajectory refinement in the design cycle of Figure 3, the uncertainty 210 in the 
models at the specific conditions may be reduced to the limit of available resources 
at design time. 

Once the uncertainty 210 is reduced during design time 220 and 
manufacturing time 230, the information gathered from these processes is used to 
identify patterns of uncertainty in each individual plant. These patterns are then 
used to sub-select model controllers at manufacturing time 230 and to further _ 
reduce uncertainty during run time 240 of each plant individually. 
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At design time 220, the highest amount of uncertainty exists about a system. 
However, the greatest amount of resources for reducing the uncertainty is also 
available at design time 220. For example, at design time 220, hours of computing 
time and large workstations are available to generate control models. The amount 
of uncertainty at manufacturing time 230 decreases, but so does the resources for 
reducing the uncertainty. For example, resources available at manufachiring time 
230 are minutes of compute time and smaller microprocessors that characterize the 
drive and find the closest matching control models from those generated at design 
time 220. 

At run time 240 the uncertainties are reduced to a lower level (compared to 
the uncertainty at design time). But the resources avaUable to reduce or handle the 
uncertainty at run time 240 are also lessened. For example, the resources that are 
avaUable at run time 240 are on the order of mOliseconds of compute time and 
digital signal processors within an individual plant (or unit) that are used to 
implement a model controller. As such, more calculations are performed at tlrnBs 
when more compute power is available. A metric of uncertainty (given the level of 
available information) is maintained and used in selecting an appropriate control 
model with a specific robustness-performance tradeoff. 

Figure 3A illustrates an embodiment of a design cycle. In one embodiment, 
the design cycle 300 may be used to develop a servo controller for a disk drive 
system. In step 310, the control problem is formulated for the system. In 
formulating the control problem, various operating conditions for performing 
evaluation and testing may be established. In one embodiment, the conditions (e.g., 
vibration and thermal) that are established are those conditions that the system is 
expected to perform under, referred to as comer conditions. 

At step 320, initial trajectory design is performed to identify the excitation 
trajectories to be used during data acquisition. In one embodiment, one or more 
prototype units of the system are used to perform the initial trajectory design. In 
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another embodiment, an initial simulation model may be used to perform the initial 
trajectory design. 

In one embodiment, the trajectories are time varying sinusoidal wave-forms 
starting near the rotation frequency of the drive. In another embodiment, the 
system is excited in a frequency range of approximately 10 Hz to the Nyquist rate of 
approximately 10 kHz. The resulting position signal (Ypos) excursions should be 
large enough so that the system is designed to spend only a small portion of time in 
regions of phase space where small stresses (during the manufacture of the disk 
drive) dominate the torque forces applied by the VCM on the disk drive head. 
These torque forces may be, for examples: micro-welding and/or pre-rolling 
friction stresses. The input signal and Ypos variations created by test conditions 
after manufacturing are also considered. These test conditions may correspond to a 
variety of mechanical and environmental tests that a manufactured drive may be 
subjected to in actual use. 

At step 330, closed loop data is acquired on the system using a prototype 
unit, as illustrated in Figure 3B. A specific operating condition is selected under 
which closed loop data acquisition is to be performed on the prototype unit. The 
trajectory 333 from step 320 of Figure 3 A, and the feedback from controller 332 is 
input to plant 331 while the corresponding control signal 334 and Ypos signal 336 
data is collected. In one embodiment, controller 332 may be a baseline controUer. 
In another embodiment, controller 332 may be a controller selected from a previous 
design iteration. 

The data from step 330 is then used to perform system identification, at step 
340. The identification generates a control model of the prototype unit based on 
characterization data of the behavior (i.e., input and output signals) of the 
prototype unit. The control model represents the plant model of the unit, including 
disturbances (as discussed in detail below with respect to Figure 4), and model 
uncertainty (additive, multiplicative and/or structural singular value bound used 
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for control design). In one embodiment, a software design tool may be used for 
system identification, for examples: Matlab* Control Module and System 
Identification Module, Matlab V-Analysis and Synthesis Toolbox, Matlab^System 
Identification Module (all available from Mathworks, Inc. of Natick, 
Massachusetts), and MATRDCx* Robust Control Module and System Identification 
Module (available from Integrated Systems, Inc. of Sunnyvale, California). 

At step 350, a set of algorithms is used to generate model controllers from the 
control models of step 340. In one embodiment, the algorithms may be generated 
using H-infinity control design methods. In other embodiments, other control 
design methods may be used, for example, ^-synthesis control design, as described 
m n-Analysis and Synthesis Toolbox User's Guide, Matlab, Mathworks Inc Natick 
Massachusetts, 1998. These control design methods are well known in the art; 
accordingly, a more detailed description of these methods is not provided herein. 

At step 360, a particular model controller, developed from control design 
step 350, is implemented in one or more prototype units and exposed to a test 
condition. In one embodiment, the model controller includes a combination of an 
estimator and controUer designed using a given control model (indexed by i as 
discussed in detail below with respect to Figure 6). A test condition may include 
one of a variety of operating conditions that a manufactured drive may experience 
The operation of the unit is characterized and its performance is diagnosed. 

Figure 3C iUustrates one embodiment of control design. In one embodiment 
a set of linearized differential equations is computed based on data from 
identification 340 of Figure 3A. At step 351, a high order model is developed as 
described below in relation to Figure 4. 

In one embodiment, at step 355, model reduction is performed by computing 
an approximate model of the linear high order model, of step 351, with a fewer 
number of states than the high order model. This may also provide an estimate of. 
errors due to the approximation. At step 356, estimator controller or controller 
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synthesis may be performed including computing an estimator and /or controller 
that optimizes certain properties, for examples, H,, Hoo, or the n norm of the closed 
loop system. At step 357, the high order model is evaluated by computing 
multivariable frequency responses, step responses, responses to stochastic noise 
processes, eigenvalues, and/or measures of robustness for the closed loop system 
consisting of the high order model and the synthesized estimator controller or 
controller. 

In another embodiment, at step 352, full order estimator controller or 
controller synthesis may be performed by computing an estimator controller that 
optimizes the properties of the closed loop system. At step 353, controller reduction 
is performed by computing an approximate estimator controller of the full order 
estimator controller of step 352, with fewer states than the full order estimator 
controller. At step 626, the high order model is evaluated by computing 
multivariable frequency responses, step responses, responses to stochastic noise 
processes, eigenvalues, and/or measures of robustness for the closed loop system 
consisting of the high order model and the estimator controller or controller from 
step 353. 

In one embodiment, at step 358, a full nonlinear simulation may be 
performed by using the estimator controller or controller from step 353 or step 356 
and running non-linear closed loop simulations for various initial data values and 
perturbed model parameters. The resulting estimator controller or controller 
satisfying robustness and performance criteria is the model controller for the 
control model. These model controllers may also be referred to as robust- 
performance design based model controllers. 

Referring back to Figure 3A, in one embodiment, the steps of design cycle 
300 may be executed multiple times to perform higher level iterations involving 
refined identification excitation trajectory, identification, and control design _ 
implementation, step 370. In one embodiment, computation of bounds on model 
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errors for identification and refinement of identification excitation trajectory may be 
performed as is disclosed in co-pending U.S. patent appUcation entitled "MODEL 
ERROR BOUNDS FOR IDENnHCATION OF STOCHASTIC MODELS FOR - 

CONTROL DESIGN," AppUcation No. , filed on , assigned to the assignee 

of the appUcation herein, and incorporated herein by reference. 

In one embodiment, tiie same prototype unit is used to repeat steps 330 to 
360 under a different operating condition. In another embodiment, additional 
prototype units may be used to repeat steps 330 to 360 in order to obtain a 
representative sample of in-use variations as well as a smaUer sample of 
manufacturing variations. In an alternative embodiment, steps 350 and 360 may 
not be performed if step 340 provides a sufficiently accurate identification of the 
system and a mathematically compatible bound on model uncertainty, for example, 
as disclosed in: co-pending U.S. patent appUcation entitied, "REAL-TIME 

PLANNER FOR DESIGN," AppUcation No. , filed on , assigned to ti.e 

assignee of tiie appUcation herein, and incorporated herein by reference; and co- 
pending U.S. patent application entitled "MODEL ERROR BOUNDS FOR 
IDENTIHCATION OF STOCHASTIC MODELS FOR CONTROL DESIGN," 
Application No. , filed on . 

In an alternative embodiment, design cycle 300 may choose different 
methods for design of a control model depending on specific design objectives as 
described, for example, in U.S. Patent No. 5,880,959 of Shah, et al., , and "REAL- 
TIME PLANNER FOR DESIGN," Application No. , filed on , both assigned 

to the assignee of the application herein and botii incorporated herein by reference. 

Figure 4 iUustrates an embodiment of a control model structure for system 
identification performed at design time. Available computing resources at design 
time 220 of Figure 2 are taken advantage of by pre-computing a large set of control 
models for different plant conditions that can occur in-use and in the _ 
manufacturing of the system. Identification with specific model structiire, as shown 
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in Figure 4, is used to obtain weU-tuned control models from limited closed loop 
data. Identification for plant model 400 allows for the affects of the three 
components: controllable observable modes 410, uncontroUable observable modes 
420, and quasi-static oscillations 430. Concurrent estimation of all components in 
the control model allows accurate estimation of each component. 

In one embodiment, identification for a servo control system may be 
performed by running closed loop tests on a prototype unit, as discussed above in 
relation to Figure 3B. The plant model 400 includes controUable and observable 
modes 410 as well as uncontrollable and observable modes 420. Plant model 400 
may be excited with uncontroUable observable modes 420 due to stochastic 
disturbances 425, for example, run-out. Run-out is a condition where a head 
deviates from a track being foUowed and may be caused by factors such as 
mechanical vibrations and electrical noise in the measurement circuitry of the head. 
These disturbances permeate through the resonant structures that are inherent in a 
disk drive. In addition, plant model 400 is excited with quasi-static oscillations 430, 
for example, quasi-periodic run-out. Quasi-periodic run-out is a slowly changing 
periodic signal that may be caused by, for examples: deformation of a disk platter 
and/ or an off-axis disk. 

In one embodiment, high order methods of estimation may be used, for 
example, auto regressive with exogenous inputs (AKK) estimation and prediction 
error methods (PEM) (See Ljung, L., System Identification, Theory for the User, 2nd 
Edition, Prentice-Hall, 1999). Alternative methods of estimation of parameters are 
well known in the art; accordingly, a more detailed description of these methods is 
not provided herein. 

Figure 5 illustrates one embodiment of an uncertainty model. In one 
embodiment, the uncertainty model 510 includes all drives operating under all 
conditions in use. As previously discussed, at design time 220 of Figure 2, it is _ 
possible to obtain in-use uncertainty for prototype drives as well as for 
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representative drives. The model for a drive under a specific operating condition is 
shown in 527 along with its uncertainty model 528. The model for that drive under 
a nominal operating condition is shown as 525. Uncertainty for aU operating 
conditions for that drive is shown as 520. 

Models for each representative drive under all operating conditions and its 
corresponding uncertainty model (e.g., model 530) may also be obtained. The 
operating conditions may include, for examples: thermal, vibration environment, 
head, platter and zone on platter. The metaric chosen to represent model variations 
is defined in terms of weighted ft-equency responses, where the weights are chosen 
to reflect sensitivity of closed-loop performance to model variations. In addition, a 
lower performance, highly robust conhroller is computed for instances when in-use 
data that determines the actiial behavior of a plant is not yet available. 

At manufactiiring time 230 of Figure 2, the drives are manufachired to a 
target specification developed at design time 220 of Figure 2. A particular 
manufactured drive may then be evaluated to obtain its nominal conta-ol model. In 
one-embodiment, off-line calculations are performed on a manufactiired drive using 
available compute resources to obtain patterns of uncertainty in order to find the 
best pre-computed (at design time 220) models around the drive to represent its in- 
use variations. By characterizing the in-use model variations around the nominal 
model for each representative drive at design time 220 of Figure 2, the in-use 
variation of a manufactiired drive may be extrapolated from its nominal model 
obtained at manufactiiring time 230 of Figure 2. The uncertainty model, thus, 
allows for fuU manufactiiring time adaptation of the servo controller for each 
production unit. This may improve manufactiiring yield and unit reliability 
because of designed-in servo tolerances to manufactiiring variations. 

At manufactiiring time 230 of Figure 2, a subset of the control models may be 
selected from tiie larger set of conti-ol models, developed at design time 220 of _ 
Figure 2, based on measured unit-to-unit manufactiiring variations. A group of 
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model controUers, corresponding to the subset of control models selected at design 
time 240 to manage extrapolated in-use variations, is implemented in a drive. A 
particular model controller within the model controller group may then be selected 
at run-time based on data derived from the operation of the drive. In one 
embodiment, the model controller group is implemented in a servo system 
controller using hardware. In another embodiment, the model controller group is 
implemented in a servo system using software. 

The control scheme described herein, allows for the uncertainty of control 
systems to be managed by using successive uncertainty identification and 
reduction, knowledge of uncertainty in the control approach, and matching 
algorithms with available compute resources. The controUer architecture is multi- 
scale adaptive because it uses successive reduction of uncertainty at different time 
scales (i.e., design, manufacturing, and in-use) and matches computational 
algorithms with available compute resources. 

Figure 6 illustrates an embodiment of a multiple model controller. In one 
embodiment, adaptive parameter model controllers are developed at design time 
220 of Figure 2. In one embodiment, if RRO is expected to vary, feedforward 
adaptation for RRO may be performed at initialization of run time 240 of Figure 2. 
Typically RRO adaptation uses the input-output response given by model i 
controllers chosen by hypothesis testing, and finds the periodic components that 
best describe the offending harmonics. The feedfoward compensates for the - 
offending harmonics by inverting the plant response at the harmonics and feeding 
them to the plant input. 

At design time 220, each model controller may be designed to handle a 
narrow range of resonance and disturbance for a specific operating condition. Then 
at manufacturing time 230 of Figure 2, a subset of model controllers 620 (Model-1 
controller to model-n controller, where n is an integer), developed from their _ 
corresponding control models, may be selected. The model controller group 620 
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selected at manufacturing time 230 of Figure 2 may be able to accommodate a range 
of resonance and disturbance that is larger than the resonance and disturbance that 
may be handled by the best single model controller. 

A controller selector 630 is used to switch between model controllers 620 
during operation of the drive using data obtained from monitoring the operation of 
the system. For example, transients may occur due to seek and settling of a head, 
generating transient signals. These transient signals may be monitored to 
determine the appropriate controller to select from among model controllers 620 
with each of the controllers designed to handle a particular transient signal range. 
A sufficient set of data is collected before switching between model controllers 620 
to reduce the probability of an incorrect model controller being selected. In one 
embodiment, parametric adaptation is performed regardless of which model 
controller is selected, however, the adaptation may be based on the states of the 
estimator within the selected model controller. 

A data evaluator 640 monitors the transient in Ypos 622 and determines if 
there is enough information to select a particular model controller and, if so, which 
model controller to select. In one embodiment, the data evaluator 640 is a model 
hypothesis tester and test confidence indicator. If the test confidence indicator 
shows a high confidence level, the hypothesis tester determines which of the model 
controllers 620 should be used. For example, the control input 626 to data evaluator 
640 may be a VCM current or actuator voltage. The VCM current or actuator 
voltage must have sufficient power at frequency bands over which dish^rbances are 
to be rejected. If there is insufficient power, the hypothesis tester may not be able to 
obtain sufficient information to determine which of the model controUers 620 is 
suited to provide the best performance. 

In one embodiment, the currently selected model controller continues to be 
used in the event that the hypothesis tester has insufficient information. In anoth^ 
embodiment, if the performance of the currently selected model controller is 
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determined to be poor and there is iiisufficient data to select an appropriate model 
controller, then a fall-back model controller 650 is selected. Fall-back model 
controller 650 is designed to be robiist enough to maintain the system untU 
sufficient data is obtained to select a more appropriate model controller. 

As previously discussed, model uncertainty may be obtained for a given 
identified model at design time 220 of Figure 2. In one embodiment, the model 
uncertainty provides botmds as an additive imcertainty for both input-output 
behavior and disturbance-output behavior. The imcertainty boimds are in a form 
so that they work with standard control design methods, for examples, H-inf inity 
robust control and ^-synthesis robust control design methods. In one embodiment, 
the control models are designed using parametric adaptation methods in which the 
input parameters to the control model are gradually changed. Typical adapted 
parameters may be, for examples: amplitudes and phase harmonics in quasi- 
periodic run-out, loop gains, and known dependencies of a model controller on a 
particular data location on a disk. 

In one embodiment, the model hypothesis testing may be performed using a 
Generalized Likelihood Ratio Test (GLT) for M-ary hypothesis. Given M disjoint 
states of the world given by Hypotheses H„, H„ H„ choose hypothesis (where 
i is an integer) that maximizes the log likelihood function X(z)=log(p(z I H,)) where 
the observed data z has a conditional probability distribution function p(z I H). 

A GLT is employed where the observed data (e.g., Ypos) also depends on 
some unknown parameter x. The unknown parameter may be, for example, 
unmeasured states of the dynamic system. In this case, the GLT consists of 
choosing hypothesis H, when max, log(p(z I x, H,)) > max, log(p(z I x, H^)). For 
model controllers with unknown parameters 9 belonging to a discrete set {0„ 9^, ... 
ej, a bank of conditional Kalman Filters yielding conditional innovations 
sequences vi = v,^ are used. The power in each conditional innovations sequence-fe 
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computed and the one with the least power indicates the chosen control model. See 
also, Anderson and More, Optimal Filtering, Prentice-Hall, 1979, pages 267-273, and 
Willsky, A., Generalized Likelihood Ratio Approach to State Estimation in Linear 
Systems Subject to Abrupt Changes, Proceedings of IEEE CDC, 1974, pages 846-853. 

Figure 7 iUustrates an embodiment of a hypothesis tester and test confidence 
indicator. Identified models (e.g., ARX) at design time 220 of Figure 2, are 
converted to their state-space innovations forms using techniques well-known in 
the art, as described in Ljung, L., System Identification, Theory for the User, 2nd 
Edition, PrePtice-HaU, 1999. The model controllers are switched by model index 
635 of Figure 6. In one embodiment, the innovation form for model i is given by: 

yt=c'xi+Ei 

where y, is the measurement at time k; and 

[Uactuatork T . 
Uperiodick J 

x'^ is the state vector at time k for model i. 

{ A\B',C\K'^ } are matrices describing model i in the state-space 

innovations form, "nie computation of coefficients are computed at design time 220 
of Figure 2. 

In one embodiment, the conditional residuals are computed at run time 240 
of Figure 2. The conditional residuals computation 710 is performed using the 
innovations models. In one embodiment, the conditional residuals, v; , are 
computed using the following equation: 
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The conditional residuals are then normalized, 720, where = Q'^^v^ are 
the normalized residuals. The normalized conditional residuals are filtered using 
weighting filters W(<») 730. These weighting filters 730 are those used for the metric 
of model variations and tiiey reflect sensitivity of closed-loop performance to model 
variations. A power estimate (I) is computed with lowpass filter 740 which receives 
magnitude square filtered residuals from magnitude square operator 735. 

In one embodiment, to perform model selection, given a previously (during 
previous sample time) selected model controller k, choose a new model controller j, 
where j does not equal k, if: 

j = arg minjflj, and 

test confidence indicator = high. 

Figure 8 illustrates an embodiment of a test confidence indicator. Test 
confidence indicator 800 includes a magnitude square operator 810, a lowpass filter 
820, a AND gate 850, and comparators 830 and 840. An input (u) 805 is applied to 
magnitude square operator 810, for example, an achiator contarol input 626 of Figure 
6. The magnitude squared of input 805 is applied to lowpass filter 820. The filtered 
output 825 is appUed to comparator 830. Comparator 830 determines if filtered 
output 825 is greater than a threshold value (THRESH,). In one embodiment, the 
threshold value of comparator 840 is between 3* i?" and 10*17 Tt 

should be noted, that may be a lowpass filtered actuator input energy 

(squared value) during tirack following determined by simulations performed at 
design time 220 of Figure 2. In another embodiment, THRESH, may be other 
values. 

A monitored parameter value from a current sample time (I.) 842 and a 
previous sample time (IJ 844 is applied to comparator 840, for example, position 
sensor signal 622 of Figure 6. Comparator 840 determines if the difference between 
input 842 and 844 is greater than a threshold value multiplied by input 844. In one 
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embodiment, the threshold value (THRESH,) is between 1.1 and 2.0. In another 
embodiment, THRESH, may be more than 2.0 or less than 1.1. 

The outputs 852 and 854 of comparators 830 and 840, respectively, are 
applied to AND gate 850 with the output 860 of AND gate 850 being the confidence 
value. In one embodiment, the lowpass filter in test confidence indicator 800 and 
filter 740 of Figure 7 are chosen such that slow variations on the order of seconds 
pass-through. In an alternative embodiment, the test confidence indicator of data 
evaluator 640 of Figure 6 may be implemented by counting for a fixed number of 
track changes larger than a predetermined number of tracks. 

As previously discussed, the controllers are switched at run-time based on 
data derived from the operation of the drive, for example, taransients due to seek 
settiing and the confidence of that data. In this manner, the system may provide a 
highly tuned controller using both parametric adaptation within a model controller 
and system adaptation. The system is adaptive by enabling the conti-ol system's 
signal path to be changed to a different contiroller based on the value of the data 
derived firom in-use operation of the drive. By using an adaptive servo control 
system, the bandwidth of the servo system may be increased. A higher bandwidth 
may yield tighter head position contarol and, thereby, increase the ti-acks per inch 
(TPI) over non-adaptive servo systems. The increased bandwidth may also result in 
improved settiing time of the system. 

As previously noted, the control scheme describe herein may be 
implemented with a wide variety of systems. In an alternative embodiment, the 
contirol scheme may be used for adaptation to unmeasured variables in non- 
resonant systems, for example, changes in viscous friction in a motor drive control. 
A motor drive conb-ol may be described by the equation: 

70 + 50 = 7,^„,, where 
where J is the inertia of tfte drive system, B is the unknown viscous friction, 9 is the 
angular position of the motor shaft, and T,,_, is the applied exogenous torque 
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generated by the electromagnetic field. In one embodiment, it is assumed that the 
viscous friction changes because of the temperature of the bearings and the 
humidity of the air, ?.nd that there are no direct measurements of the bearing 
temperature or the air humidity. For a fixed motor position control system that 
uses speed measurement, changing viscous friction would deteriorate position 
control performance. The resulting system performance would be sensitive to 
immeasiured temperature and humidity variations. An adaptive system that 
estimates changing djmamics of the motor drive system and adapts the controller 
parameters to the changes would allow the control system performance to remain 
constant in spite of variatioris in the djmamics. 

In the foregoing specification, the invention has been described with 
reference to specific exemplary embodiments thereof. It will, however, be evident 
that various modifications and changes may be made thereto without departing 
from the broader spirit and scope of the invention as set forth in the appended 
claims. The specification and drawings are, accordingly, to be regarded in an 
illustrative rather than a restrictive sense. 
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CLAIMS 

What is claimed is: 

1. A method of generating a control system, comprising: 
generating a plurality of control models of the plant based on an 

operation of a prototype at a design time; 

characterizing an operation of the plant at a manufacturing time; and 
selecting a plurahty of model controllers based on the characterizing 
of the operation of the plant at manufacturing time, each of the plurality of 
model controllers corresponding to one of the plurality of control models. 

2. The method of claim 1, further comprising selecting one of the 
plurality of model controllers at a nm time for the plant based on an 
operating parameter. 

3. The method of claim 1, wherein the plurality of control models 
includes an input-output model, disturbance output stochastic model, and 
quasi-periodic disturbances. 

4. The method of claim 1, wherein generating the plurality of control 
models comprises pre-computing a plurality of model controllers based on 
robust-performance design based model controllers for different plant 
conditions that occur in-use and in-manufacturing of the plant. 

5. The method of claim 4, wherein generating the plurality of control 
models at the design Hme further comprises developing an uncertainty 
model. 

6. The method of claim 5, wherein developing an uncertainty model 
comprises: ~^ 
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obtaining a nominal model vmder a specific operating condition based 
on the operation of the prototype; and 

obtaining a plurality of model variations of the prototype. 

7. The method of claim 6, wherein developing an uncertainty model 
further comprises choosing a metric to represent the plurality of model 
variations. 

8. The method of claim 7, wherein the metric is defined in terms 
weighted frequency responses. 

9. The method of claim 8, wherein the metric reflects sensitivity of 
closed-loop performance to variation in a plant model. 



are 



10. The method of claim 5, wherein the plurality of control models 
excited with trajectories in a frequency range, each of the plurality of control 
models designed to function in a portion of the frequet^ry range. 

11. The method of claim 2, wherein characterizing the operation of the 
plant comprises measuring unit-to-unit variations and in-unit variations in 
the operation of the plant. 

12. The method of claim 11, further comprising implementing the 
pluralit>= of model controllers in a system controller. 



13. The method of claim 1 1, wherein the plurality of model controllers 
hardware implemented in the system controller. 

14. The method of claim 11, wherein the pluraUty of model controllers 
software implemented in the system controller. 

15. A controller, comprising: 
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a pluraUty of model controUers, each of the plurality of model 
controllers coupled to receive a pluraUty of monitored signals, each of the 
plurality of model controllers to generate a control signal; and 

a data evaluator coupled to receive the plurality of monitored signals 
and generate an output, the data evaluator to determine when sufficient 
information from the plurality of monitored signals is received, the data 
evaluator coupled to select the control signal of one of the pluraUty of model 
controllers based on the output of the data evaluator. 

16. TTie controller of claim 15, wherein the data evaluator comprises a 
model hypothesis tester and test confidence indicator. 

17. The controller of claim 16, wherein the test confidence indicator 
utiUzes a generaUzed likelihood ratio test to determine when sufficient 
information from the plurality of monitored signals is received. 

18. Tb.e controller of claim 16, hirther comprising a fall-back controUer 
coupled to receive the plurality of monitored signals, the fall-back controller 
generating the control signal and wherein the data evaluator is coupled to 
select the control signal of the faU-back controller when insufficient 
information from the pluraUty of monitored signals is received. 

19. The controller of claim 18, further comprising a controller selector 
coupled to receive the output of the data evaluator, the controller selector 
coupled to receive the control signal of each of the plurality of model 
controllers and the fall-back controller, the controller selector to output the 
control signal of one of the plurality of model controllers and the fall-back 
controller based on output of the data evaluator. 
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20. The controller of claim 19, wherein the plurality of model controllers 
have adaptive parameters. 

21 . The controller of claim 20, wherein the parameters of the plurality of 
model controllers are gradually adjusted. 

22. The controller of claim 21, wherein plurality of monitored signals are 
transients in a head position signal in a disk drive system and wherein the ' 
control signal adjusts a current in a motor in Hie disk drive system. 

23. A method of controlling a plant, comprising: 

receiving a plurality of monitored signals by a plurality of model 
controllers and a data evaluator; 

evaluating the plurality of monitored signals to determine when 
sufficient information of the plurality of monitored signals is received; and 

selecting one of the plurality of model controllers to generate a control 
signal based on the plurality of monitored signals. 

24. The method of claim 23, wherein each of the plurality of the model 
controllers have adaptive parameters affecting the generated control signal 
and wherein the method further comprises gradually adjusting the adaptive 
parameters of the selected one of the pluraUty of model controllers. 

25. The method of claim 23, wherein evaluating comprises performing a 
generalized likelihood ratio test on the plurality of monitored signals to 
determine which of the plurality of model controllers to select. 

26. A controller, comprising: 
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a pluraUty of model controUers, each of the plurality of model 
controllers coupled to receive a plurality of monitored signals, each of the 
plurality of model controllers to generate a control signal; and 

a data evaluator coupled to receive the plurality of monitored signals 
and generate an output, the data evaluator to determine when sufficient 
information from the plurality of monitored signals is received, the data 
evaluator coupled to select the control signal of one of the plurality of model 
controllers based on a generalized Ukelihood ratio test with frequency 
weighting to reflect closed-loop performance sensitivity to variations in a 
plant model. 

27. The controller of claim 26, wherein the data evaluator comprises a test 
confidence indicator to delay switching among the plurality of model 
controllers imtil confidence in a hjrpothesis test is high. 

28. Hie controller of claim 27, wherein the confidence is based on a 
number of seek settles in a given time. 

29. The controller of claim 27, wherein the confidence is based on an 
average input power being high. 

30. The controller of claim 27, wherein the confidence is based on a 
difference in log-likelihood ratio estimates. 

31. The controller of claim 27, further comprising a fall-back controller 
coupled to receive the pluraUty of monitored signals, the fall-back controller 
generating the control signal and wherein the data evaluator is coupled to 
select the control signal of the fall-back controller when confidence in the 
hypothesis test is low. 
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32. A disk drive system, comprising: 
a head assembly; and 

a servo controller coupled to the head assembly, the servo controller 
comprising: 

a plurality of model controllers, each of the plurality of model 
controllers coupled to receive a plurality of head position sensor 
signals, each of the plurality of model controllers to generate a control 
signal; and 

a model hypothesis tester and test confidence indicator coupled 
to receive the plurality of head position sensor signals and generate a 
model index, the model hypothesis tester and test confidence 
indicator to determine when sufficient information from the plurality 
of head position sensor signals is received, the model hypothesis 
tester and test confidence indicator coupled to select the control signal 
of one of the plurality ef iaodel controllers based on the model index. 

33. The disk drive system of claim 32, wherein the servo controller further 
comprises: 

a fall-back controller coupled to receive the plurality of head position 
sensor signals, the fall-back controller generating the control signal and 
wherein the model hypothesis tester and test confidence indicator is coupled 
to select the control signal of the fall-back controller when insufficient 
information from the plurality of head position sensor signals is received; 
and 

a controller selector coupled to receive the model index, the controller 
selector coupled to receive the control signal of each of the plurality of model 
controllers and the fall-back controller, the controller selector to output th^ 
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control signal of one of the plurality of model controllers based on the model 
index. 

34. The disk drive system of claim 33, wherein the control signal is a voice 
coil motor current. 

35. The disk drive system of claim 33, wherein the control signal is an 
actuator voltage. 

36. A method of handling disk drive variatioi\s, comprising: 
receiving a pluraUty of position sensor signals and a control signal by 

a servo controller, the plurality of position sensor signals having a range of a 
factor; 

evaluating the plurality of position sensor signals to determine when 
sufficient information of the plurality of position sensor signals is received; 
and 

selecting one of a plurality of model controllers within the servo 
controUer to operate on the control signal based on the factor. 

37. The method of claim 36, wherein evaluating comprises performing a 
generalized likelihood ratio test on the plurality of position sensor signals to 
determine when there is a high confidence in the plurality of position sensor 
signals received by the servo controller. 

38. The method of claim 37, wherein each of the plurality of the model 
controllers have adaptive parameters affecting the received control signal 
and wherein the method further comprises gradually adjusting the adaptive 
parameters of the selected one of the plurality of model controllers. 

39. The method of claim 36, wherein the factor is resonance. — 
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40. The method of claim 36, wherein the factor is disturbance. 

41. The method of claim 40, wherein the factor is environment. 

42. A method of designing control models, comprising: 
manufacturing a first prototype of a system; 

characterizing the first prototype of the system to obtain a first data; 

and 

computing a plurality of control models based on the first data of the 
system. 

43. The method of claim 42, further comprising: 
generating a computational simulation model of the system; 
operating the simulation model to obtain a second data of the system; 

and 

computing the plurality of control models based on the first data and 
the second data of the system. 

44. The method of claim 43, further comprising: 
manufacturing a plurality of prototypes of the system; 
characterizing manufacturing variations of the plurality of prototypes 

of the system to obtain a third data; and 

computing the plurality of control models based on the first data, the 
second data, and the third data. 

45. The method of claim 43, further comprising performing a plurality of 
iterations of refining the simulation model and computing the plurality of 
control models. 

-31- 



wo 00/41042 



PCT/USOO/00387 



46. The method of claim 45, wherein computing the plurality of control 
models comprises exciting the plurality of control models with trajectories at 
an operating condition, each of the plurality of control models designed to 
function at the operating condition. 

47. A method of manufacturing disk drives, comprising: 

designing a plurality of control models and a corresponding plurality 
of uncertainty models of a servo system imder at least one operating 
condition; 

manufacturing a disk drive; and 

characterizing the disk drive to obtain a nominal control model from 
the plurality of control models. 

48. The method of claim 47, further comprising implementing a pluraUty 
of model controllers corresponding to a set of the plurality of control models 

Jr. a -controller of the servo system based on the nominal control model. 

49. The method of claim 48, further comprising selecting one of the 
plurality of model controllers during operation of the disk drive. 

50. The method of claim 48, wherein the at least one operating condition 
is a condition selected from a group consisting of temperature, humidity, 
vibration spectrum, and age. 
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date and not in conflict with the application but cited to undersund 
the principle or theory underlying the invention 

document of particular relevance; the claimed invention cannot be 
considered novel or cannot be considered to involve an inventive step 
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3- Claims Nos.: 
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1. [x] ^^-^-<«^tiona,sea„^feesweretimelypaidbytheappUcam.U^ 

° ^Z'^T^ IT '^'^^ effort justifying an additional fee. this Authority did 
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4- No required additional search fees were timely paid by 



restricted to the invention first mentioned in the cl 
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[x] The additional search fees were accompanied by the applicant' s protest 
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BOX II. OBSERVATIONS WHERE UNITY OF INVENTION WAS LACKING 
This ISA found multiple inventions as follows: 

This application conlains the following inventions or groups of inventioDs which are not so linked as to form a 
single inventive concept under PCT Rule 13. K In order for all inventions to be searched, the appropriate additional 
search fees must be paid. 

Group I, claim(s) 1-31 and 42-46. drawn to a controller, a control system and a method of designing, 
generating conirol system for controlling a plant. 

Group II. claim(s) 32-41, drawn to disk drive system and a method of handling a disk drive variation, . 
Group III, claim(s) 47-50, drawn to a method of manufacturing disk drives. 

The inventions listed as Groups I, II and III do not relate to a single inventive concept under PCT Rule 13. 1 
because, under PCT Rule 13.2, they lack the same or corresponding special technical features for the following reasons: 
group I relates to control system and a method of controlling a control system being examined in US class/subclass 
700/19. and Group II relates to a disk drive system and a method of handling a disk drive variation being examined in 
class/subclass 360/75 and Group III relates to method of manufaauring disk drives being examined in class/subclass 
700/117. All three grc*^ps as stated above are being examined by different class and/ or subclass and they do not relate 
to a single inventive concept under CT Rule 13.1. 
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